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This  report  describes  a  method  for  using  three-dimensional  infinite  elements  to  compute  acoustic  fields  in  the  infinite 
region  containing  water  that  surrounds  a  submerged  structure.  In  this  method,  finite  elements  are  used  to  represent  a 
bounded  region  that  contains  the  structure  a„d  ir  some  cases  contains  part  of  the  surrounding  water.  The  unbounded  region 
that  surrounds  the  bounded  region  represented  by  finite  elements  and  that  contains  the  remainder  of  the  water  is  represented 
with  infinite  elements.  The  infinite  elements  implement  the  added  mass  approximation  to  the  Helmholtz  equation. 

Several  other  methods  are  available  for  modeling  the  effects  of  the  surrounding  water.  For  the  added  mass  approxima 
tion,  an  alternative  io  using  infinite  elements  is  to  model  parts  of  the  suirounding  water  with  several  layers  of  finite  ele¬ 
ments.  To  determine  whether  the  method  using  infinite  elements  gives  an  advantage,  the  effieieney  of  that  method  is  com 
pared  with  the  efficiency  of  the  alternative  method.  The  relative  efficiencies  ,..e  determined  by  comparing,  for  each  method, 
the  effort  to  construct  the  model  and  the  accuracy  of  the  solutions  obtained.  Evaluation  of  the  infinite  elements  shows  that 
they  produce  accurate  results  especially  if  they  are  used  over  a  layer  of  finite  fluid  elements.  However,  when  solutions— 
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ABSTRACT 


This  report  describes  a  method  for  using  three-dimensional 
infinite  chments  to  compute  acoustic  fields  in  the  infinite  region 
containing  water  that  surrounds  a  submerged  structure.  In  this 
method,  finite  elements  are  used  to  represent  a  bounded  region  that 
contains  the  structure  and  in  some  cases  contains  part  of  the 
surrounding  water.  The  unbounded  region  that  surrounds  the 
bounded  region  represented  by  finite  elements  and  that  contains  the 
remainder  of  the  water  is  represented  with  infinite  elements.  The 
infinite  elements  implement  the  added  mass  approximation  to  the 
Helmholtz  equation. 

Several  other  methods  are  available  for  modeling  the  the  effects 
of  the  surrounding  water.  For  the  added  mass  approximation,  an 
alternative  to  using  infinite  elements  is  to  model  part  of  the 
surrounding  water  with  several  layers  of  finite  elements.  To 
determine  whether  the  method  using  infinite  elements  gives  an 
advantage,  the  efficiency  of  that  method  is  compared  with  the 
efficiency  of  the  alternative  method.  The  relative  efficiencies  are 
determined  by  comparing,  for  each  method,  the  effort  to  construct  the 
model  and  the  accuracy  of  the  solutions  obtained.  Evaluation  of  the 
infinite  elements  shows  that  they  produce  accurate  results  especially  if 
they  are  used  over  a  layer  of  finite  fluid  elements.  However,  when 
solutions  obtained  using  infinite  elements  are  compared  with  solutions 
obtained  using  only  finite  elements,  it  is  seen  that,  for  added  mass 
problems,  infinite  elements  do  not  offer  a  significant  advantage  over 
finite  elements.  Infinite  elements  arc  expected  to  provide  a  greater 
advantage  when  used  to  solve  the  Helmholtz  equation. 


ADMINISTRATIVE  INFORMATION 

Tins  project  was  jointly  supported  by  the  DTNSRDC  Independent  Research  Program,  sponsored  by  the 
Space  and  Naval  Warfare  Systems  Command,  Director  of  Navy  Laboratories,  SPA  WARS  005,  and 
administered  by  the  Director  of  Research,  DTNSRDC  012.3  under  The  Office  of  Chief  of  Naval  Research 
Program  Element  61 153N  (Task  Area  ROOO-Ol-Ol  under  DTNSRDC  Work  Unit  18-14-125.)  and  by  the  Office 
of  Naval  Research  Mathematical  Sciences  program,  Numerical  Methods  for  Naval  Vehicles  Program,  Program 
Element  61 153N,  Task  Area  BR-014-03-51,  Work  Unit  1840-040. 

INTRODUCTION 

To  compute  magnetic  or  acoustic  fields  about  submerged  marine  structures,  the  structure,  the 
surrounding  water,  and  interactions  between  them  must  be  modeled.  The  field  and  the  structure  can  each  be 
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represented  by  finite  elements  and  the  representations  can  be  coupled  using  a  method  developed  by 
Zienkiewicz  and  Newton.*  If  the  structure  is  submerged  in  the  sea,  the  extent  of  the  water  surrounding  it  is 
so  great  that  the  region  containing  the  water  is  best  represented  as  having  infinite  extent.  If  finite  elements 
are  used  to  solve  such  a  problem,  one  can  model  the  structure  and  a  finite  part  of  the  infinite  region 
surrounding  the  structure.  This  leaves  the  problem  of  accounting  for  the  remainder  of  the  surrounding 
infinite  region. 

One  method  of  accounting  for  the  infinite  region  of  water  surrounding  a  submerged  structure  is  to 
represent  the  region  with  infinite  elements.  To  solve  problems  involving  sound  radiated  or  scattered  by  a 
submerged  structure,  the  water  is  treated  as  an  acoustic  fluid,  that  is,  a  fluid  having  an  acoustic  pressure  field 
described  by  the  wave  equation.  For  low  frequencies,  an  approximation  called  the  added  mass  approximation 
(described  later  in  this  report)  simplifies  the  solution  of  the  acoustic  problems.  This  report  documents  the 
development  of  infinite  elements  and  the  evaluation  of  these  elements  for  fluid-structure  interaction  probler. 
involving  the  vibrations  of  submerged  structures  in  which  the  added  mass  approximation  is  made. 

To  be  useful,  an  infinite  element  must  not  only  be  able  to  produce  solutions  with  good  accuracy  but 
must  achieve  this  accuracy  while  requiring  less  effort  for  modeling  and  computation  than  competing  methods. 
Several  competing  methods  of  accounting  for  the  surrounding  infinite  region  are  as  follows: 

•  The  surrounding  region  is  truncated  by  modeling  a  large  part  of  it  with  finite  elements  and 
applying  boundary  conditions  that  approximate  the  effects  of  the  remainder  of  the  infinite 
region. 

•  The  solution  is  expanded  as  a  series  of  analytic  functions  in  the  surrounding  region  with  the 
coefficients  of  the  series  terms  introduced  as  unknowns. 

•  The  effects  of  the  surrounding  region  are  represented  by  an  integral  equation  on  the  boundary 
of  the  region  surrounding  the  finite  elements. 

O 

Zienkiewicz  et  al."  compare  modeling  with  these  methods  and  with  infinite  elements  and  find  that  each  has 
certain  advantages  and  disadvantages. 
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The  advantages  of  using  infinite  elements  are  that  they  maintain  symmetry  and  reasonable  bandwidths 

in  the  matrices  produced,  and  the  method  is  easy  to  implement  using  available  finite  element  programs.  On 

the  other  hand,  some  formulations  of  infinite  elements  introduce  more  degrees  of  freedom  than  do  competing 

methods,  and  in  some  cases  an  extra  effort  must  be  made  to  ensure  that  element  shapes  arid  modeling 

3 

configurations  provide  unique  mappings.  See  Bettess  and  Bettess  for  a  discussion  on  avoiding  nonunique 
mappings.  For  some  problems,  these  characteristics  may  require  additional  efTort  in  producing  the  numerical 
model. 


In  this  report  the  surrounding  infinite  medium  is  assumed  to  be  an  acoustic  fluid  and  infinite  elements 
are  used  to  represent  th»o  fluid.  An  infinite  element  represents  a  sector  of  the  infinite  region  extending  from 
the  boundary  of  die  finite  region.  In  both  finite  and  infinite  elements  the  unknown  function  is  approximated 
by  shape  functions  and  afunctional  involving  the  shape  functions  is  integrated  over  the  area  or  volume  of  the 
element.  To  obtain  convergent  integrals  for  infinite  elements,  one  of  two  schemes  is  used.  The  first 
incorporates  decay  factors  in  the  shape  functions  that  vary  in  the  direction  that  extends  to  infinity,  the  second 
-.cheine  maps  the  infinite  element  into  a  standard  square  or  cube.  In  the  second  scheme,  the  decay  factors 
appear  m  the  mappings  that  compress  the  infinite  element  into  a  finite  region.  In  either  case,  the  decay 
factors  may  iuke  one  of  several  forms,  they  may  decrease  exponentially  with  the  distance  r  from  a  fixed  point, 
ot  decrease  as  a  povvei  of  r.  Many  earlier  developments  of  infinite  elements  used  exponential  decay  factors,  * 
but  tiles j  require  the  choice  of  a  somewhat  arbitrary  parameter  called  the  decay  length,  and  at  moderate  radii 
the  shape  functions  do  not  approximate  the  behavior  of  solutions  as  well  as  the  reciprocal  of  a  power  of  the 
radius.  Later  work  in  infinite  elements’  has  emphasized  decay  factors  that  decrease  as  l/r",  and  these  factors 
are  used  in  this  work. 


Several  beginnings  of  applications  of  infinite  elements  to  vibiations  of  submerged  structures  have  been 


made  in  which  promising  lesults  have  been  reported.  ’  The  infinite  elements  used  in  these  papers  contained 


exponential  decay  factois,  and  the  papcis  do  not  report  useful  comparisons  with  alt*  rnative  methods  of 


computation. 
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In  problems  with  steady  stau  solutions  for  the  acoustic  response  of  submerged  structures,  the  behavior 
of  the  surrounding  acoustic  fluid  is  described  by  the  Helmholtz  equation.  For  low  frequencies,  an 
approximation  can  be  made  in  which  the  effect  of  the  surrounding  fluid  becomes  that  of  increasing  the  mass 
of  the  structure.  This  approximation  is  called  the  added  mass  approximation.  The  results  given  in  this  report 
show  that,  for  fluid-st.ucture  interaction  problems  using  the  added  mass  approximation,  infinite  elemi  its  give 
good  results  and  ,-n  fact  are  competitive  with  an  alternative  method.  Howeve  m  the  case  of  the  added  mass 
approximation  the  infinite  elements  do  not  offer  a  clear  advantage  overt'  *  alternative.  This  conclusion  does 
not  apply  to  the  application  of  infinite  elements  to  Helmholtz  equation  problems  in  which  infinr  •  elements  ar“ 
expected  to  offer  advantages  over  their  competitors,  since  for  these  problems  the  _.oi  lie  wave  characteristics 
of  the  sound  field  are  important  and  the  infinite  elements  will  reflect  these  charaeterist’cs. 


FINITE  AND  INFINITE  ELEMENT  MODEL  OF  A  FLUID-STRUCTURE  SYSTEM 
A  finite  element  model  is  used  to  compute  the  natural  frequencies  'f  a  structure  sjbmerged  in  an 
acoustic  fluid  or  the  structure's  response  to  harmonic  excitation  over  a  range  of  frequencies  The  finite 
element  model  of  the  acoustic  fluid  in  which  a  structure  is  submerged  me;  contain  both  finite  and  infinite 
fluid  elements.  Both  finite  and  infinite  fluid  elements  may  be  of  various  orders  with  various  numbers  of 
nodes.  The  infinite  and  finite  fluid  elements  used  in  this  report  are  brick-type  three-dimensional  elements 
with  eight  nodes. 

Paradoxically,  an  infinite  element  is  a  type  of  finite  element.  (The  adjective  “finite”  indicates  that  a 
finite  element  has  at  least  finite  dimensions,  that  it  is  not  an  infinitesimal  element  as  used  in  the  derivation  of 
equations  of  physics.)  Thus  infinite  elements  inherit  many  of  the  desirable  features  of  finite  elements  The 
configuration  of  the  model  in  terms  of  nodes  and  connections  and  the  form  of  the  resulting  matiices  do  not 
change  when  infinite  elements  are  included.  Computing  the  matrix  entries  contributed  by  infinite  elements 
follows  the  same  steps  as  computing  matrix  entries  for  finite  elements,  although  some  special  considerations 
are  needed  to  ensure  that  integrals  over  infinite  regions  converge.  'Die  contribution  of  an  infinite  element  is 
included  in  a  finite  element  matrix  in  the  same  manner  as  the  contribution  of  a  finite  element  The  finite 
element  matrices  remain  symmetric  and  banded  if  infinite  elements  are  included.  These  properties  make  it 
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easy  to  include  mholte  element*  in  existing  finite  element  codes,  and  in  the  following  discussions  we  will  use 
expressions  such  «s  finite  element  method,  finite  element  model,  or  finite  element  matrix,  even  though 
infin't;  elements  may  be  included. 

To  u.  e  infinite  elements  to  model  a  structure  submerged  in  an  unuounded  acoustic  fluid  such  as  the 
sea,  one  model;  x  finite  legion  containing  the  structure  and  optionally  part  of  the  fluid  close  to  the  structure 
with  finite  i  lenient*  and  the  remainder  with  infinite  elements.  A  model  of  tilt  surrounding  fluid  that  contains 
only  Unit,  fluid  elements  lunsists  of  one  or  more  layers  of  finite  elements  covering  the  surface  formed  by  the 
interface  l>  tween  the  structure  and  the  fluid.  If  infinite  fluid  elements  are  used,  one  or  more  outer  layers  of 
finite  ekmeius  are  replaced  by  a  layer  of  infinite  elements.  Each  infinite  element  represents  a  sector  of  the 
infinite  region  radiating  outward  from  the  boundary  of  the  region  modeled  with  finite  elements.  Figure  1 
ohoivo  the  three  configui  ttion*  of  finite  and  infinite  elements  modeling  the  sairounding  fluid.  The  most 
advantageous  ariungement  in  terms  of  reducing  modeling  and  computation  cosis  is  that  all  the  finite  fluid 
elements  are  replaced  so  that  only  infinite  elements  remain  to  model  the  acoustic  fluid. 

'flic  finite  elemtnt  1.  othod  prod  ices  one  matrix  equation  for  the  displacements  in  the  structure  and 
another  matrb  equation  for  the  acoustic  pressure  field  in  the  surrounding  fluid.  These  matrix  equations  are 
u  upled  by  terms  arising  from  interactions  between  the  fluid  and  the  structure.  To  compute  natural 
frequences  of  the  fluid-structure  system,  one  assumes  a  time  harmonic  solution,  and  the  matrix  equation 
Lee  amis  the  >  quation  for  an  eigenvalue  problem  that  detern  ines  the  frequencies.  To  compute  the  response 
.a  a  harmonic  exciting  lorce,  one  solves  a  system  of  equations  with  complex  coefficients  and  a  sinusoidal 
forcing  term  for  the  nodal  displacements. 

Standard  finite  elemtnt  modeling  procedures  produce  the  matrix  equation  for  the  vector  of  nodal 
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displacements  u  in  the  structure 
Mu  f  Bu  -r  Ku  —  f 

in  which  M,  B,  and  K  are  the  structural  mas*,  damping,  and  stiffness  mat.ricies  and  f  is  the  vector  of  forces 
acting  on  the  nodes  of  the  structure. 
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Fig.  1.  Infinite  elements  and  finite  elements  modeling  the  surrounding  fluid: 

(a)  Fluid  modeled  with  three  layers  of  finite  elements,  (b)  Fluid  modeled 
with  infinite  elements  over  two  layers  of  finite  elements,  and  (c)  Fluid 
modeled  with  only  infinite  elements. 


The  differential  equation  for  the  acoustic  pressure  P  in  the  fluid  is  the  wave  equation 


1  92P 


^T-^P=° 


where  c  is  the  speed  of  sound  in  the  fluid.  A  method  of  analogies  that  uses  solid  three-dimensional  structural 


elements  produces  a  matrix  equation  for  the  pressure  field  in  the  acoustic  fluid.  The  method  of  analogies  is 


applied  by  giving  Young’s  modulus  15  =  1.0  and  the  shear  modulus  0  =  lO5  and  produces  a  matrix  equation 
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for  the  pressure  field  in  the  acoustic  fluid 

Qp  +  Hp  =  0 

where  p  is  a  vector  of  pressures  at  the  nodes  of  the  fluid  finite  elements. 

If  the  acoustic  pressure  is  harmonic  in  time  with  frequency  w,  that  is,  P(x,t) 
point  in  space  and  k  =  oa/c,  the  wave  equation  becomes  the  Helmholtz  equation 

k2p  +  V2p  --=  0. 

Interactions  at  the  fluid-structure  interface  couple  the  structure  and  fluid  matrix  equations  See 


p(x)e‘“l  where  x  is  a 


(1) 


Schroeder  and  Marcus^  and  Everstine  et  al.^  for  details  of  the  analysis  of  the  fluid-structure  coupling  The 
acoustic  pressure  acting  on  the  interface  produces  a  force  vector  (with  the  convention  that  forces  acting 
outward  from  the  surface  are  positive),  so  the  matrix  equation  for  the  structure  bee  >mcs 

Mu  +  Bu  +  Ku  =  -  Ap  +  F 

where  the  matrices  M,  B,  and  K  are  defined  previously,  A  is  a  matrix  whose  entries  reflect  the  areas  of 
elements  on  the  interface,  and  F  is  the  vector  of  forces  other  than  the  acoustic  forces  on  the  nodes  of  the 
structure.  The  mo,  >n  of  the  structure  also  affects  the  acoustic  pressure  field  in  the  fluid.  The  normal 
component  us  of  the  acceleration  of  the  interface  between  the  structure  and  fluid  is  related  to  the  gradient  of 
the  acoustic  pressure  field  by  the  equation 
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where  p  is  the  density  of  the  fluid  and  d/Ov  is  the  derivative  in  the  direction  ii  normal  to  the  interface  This 
relation  is  added  as  a  boundary  condition  and  produces  a  force-type  term  in  the  equation  for  the  acoustic 
pressure  field.  In  the  fluid  finite  element  matrix  equation,  this  term  is  an  area  matrix  term 


Qp  +  Hp 


pA  r. 
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For  the  time  harmonic  case  this  equation  becomes 


k2  Qp  +  Hp  =  -  pAT.  (2) 

The  matrix  A  is  the  same  area  matrix  that  appears  in  the  matrix  equation  for  the  structure.  Either  a 
consistent  or  a  lumped  formulation  may  be  chosen  to  compute  the  matrix  A.  For  problems  using  the  added 
mass  approximation,  the  lumped  formulation  has  worked  well^’^  and  was  used  in  this  work. 

Combining  the  matrix  equations  for  the  structure  and  the  acoustic  field  produces  a  system  of  equations 
for  the  coupled  fluid-structure  problem 

Mu  +  Bu  +  Ku  +  Ap  =  F 

T-  •  (3) 

-pA^+Qp  +  Hp  =  0. 


These  equations  are  written  in  the  form  of  a  matrix  equation 
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This  equation  can  be  made  symmetric  by  defining  a  new  variable,  essentially  a  velocity  potential 
t 

q(t)  =  /  p(r)dr, 
o 

and  integrating  the  second  row  to  obtain  the  matrix  equation  in  'he  form 


M  O' 
0  Q 


+ 


'  B  A 

—  Af  0 


k  o 
+  [°  H 


If  the  equation  is  harmonic  in  time  (this  is  the  situation  to  which  the  Helmholtz  equation  applies),  then  the 
coupled  matrix  equation  becomes 


The  response  to  a  harmonic  excitation  of  frequency  w  is  computed  by  solving  the  preceding  equation  for 
displacements  at  selected  grid  points  resulting  from  an  applied  harmonic  load  F(t)  =  e'ut  applied  at  a  load 
point.  For  natural  frequencies,  the  eigenvalue  equation 

[mo]  [  b  a]  [k  o’ 

w  [0  Qj  +w[-At  oj  +  [0  Hj=  0 

is  solved  for  the  complex  frequencies  oj. 

In  this  report  the  added  mass  approximation,  valid  for  low  frequencies,  is  used.  For  small  oj  the 
Helmholtz  equation  (Eq.  l)  is  approximated  by  y2P  —  0,  the  fluid  matrix  equation  (Eq.  2)  becomes 
lip  =  -pAT;  and  the  coupled  fluid-structure  system  of  equations  (Eq.  3)  is 

Mu  +  Bu  +  Ku  +  Ap  =  0 
-pATii  +  Hp  =  0. 

If  the  second  equation  is  solved  for  p  and  this  expression  for  p  is  substituted  into  the  first  equation,  the  first 
equation  Mil  include  the  term  pAH  'At  if,  which  has  the  form  of  a  mass  term  and  is  called  the  added  mass. 
This  substitution  was  used  to  compute  the  natural  frequencies  of  a  submerged  plate.  For  computational 
reasons,  the  substitution  was  not  used  explicitly  for  the  frequency  response  of  a  submerged  cylinder. 

INFINITE  ELEMENT  FORMULATION 

5 

The  infinite  elements  used  in  this  study  are  of  the  mapped  reciprocal  decay  type.  They  are  three- 
dimensional  brick-type  elements  with  eight  nodes.  The  formulation  of  these  elements  is  similar  to  that  of  an 
isoparametric  finite  element  in  that  the  region  represented  by  the  clement  is  mapped  into  a  standard  finite 
cute.  The  mapping  for  an  infinite  element  incorporates  decay  factors,  factors  that  decrease  to  zero  fast 
enough  to  map  the  infinite  region  into  a  standard  finite  cube.  The  decay  factors  used  in  these  elements  are 
reciprocals  of  powers  of  the  distance  r  from  a  convenient  origin  to  the  mapped  point. 
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The  mapping  function  is 

x  =  x1F,($)F2(,)H,({)  +xiF1(?)F2(,)H1(s) 

+  x»Fi(«JFl(l)Hi(f)  +  x<F2(e)F,U)H1(f) 

+  xtfteWMMs)  +  XePl(«)F2(>))H2(c) 

+  x,F,(e)F1(i,)II2(f)+x8F!({)F,(,)H2(S) 

and  the  shape  functions  are 

N,(«,o,f)  =  F2({)F2(,)G,(t) 

N2(?,i|,s)  =  F,(5)F2(i,)G1(s) 

Ns(S,w)  =  F,(?)FI(,)G,(f) 

M,(S,>1,s)  =  Fj(«)F1(l|)G1($) 

Ns(?,>),{)  =  F2U)F2(,)G2(t) 

NeU.x.s)  =  FrfOFrfxJG^f) 

N7(f,»,j)  =  F1(f)F1(„)G2(t) 

N8(e,1,s)  -  F2(«)F1(7,)G2(f) 

where  the  functions  F,  G,  and  II  are  defined  by 

Pi(r)  =  (1-  r)/2;  F2(t)  =  (1+t)/2; 

Gi(?)  =  ?(?-i)/2;  G2(?)  =  i-?2; 

H,(?)  =  -2?/(i-f);  H2(?)  =  (i+?)/(i-?)- 

In  the  preceding  definitions,  a  point  in  space  is  x  =  xi  +  yj  +  zk,  for  i,  j,  and  k  unit  vectors  in  the  x-,  y-,  and 
z-directions,  and  x;  and  N;  are  the  locations  and  shape  functions  for  the  ith  node  in  an  infinite  element.  Using 
these  functions,  one  can  compute  the  Jacobian  of  the  mapping  function  and  the  gradients  of  the  shape 


functions. 


x.  ■ 
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the  infinite  element  correspond  to  the  four  nodes  in  the  standard  cube  at  which  f  —  0.  Since  the  infinite 
element  is  mapped  into  a  finite  region,  a  standard  Gaussian  quadrature  is  used  to  compute  the  preceding 
integrals. 

EVALUATION  OF  THE  INFINITE  ELEMENTS 

The  infinite  elements  developed  for  use  with  the  added  mass  approximation  .'ere  evaluated  by 

performance  tests  in  which  the  accuracy  of  solutions  obtained  using  infinite  elements  was  compared  with  the 

accuracy  obtained  using  only  finite  elements.  These  tests  of  the  infinite  elements  are  intended  not  only  to 

show  whether  the  elements  produce  accurate  solutions,  but  also  to  determine  whether  they  offer  an  advantage 

over  other  methods,  that  is,  whether  they  produce  as  good  or  better  accuracy  with  less  modeling  effort  or  less 

computer  time.  Although  both  finite  and  infinite  elements  can  be  formulated  with  higher  order  mapping  and 

shape  functions,  both  were  chosen  to  have  eight  nodes  and  to  have  shape  functions  that  are  linear  in  the  plane 

of  the  fluid-structure  interface.  The  infinite  elements  with  eight  nodes  have  shape  functions  that  are  quadratic 

in  the  infinite  direction.  For  fluid-structure  interaction  problems  using  the  added  mass  approximation,  it  has 

been  shown  that  jood  results  can  be  obtained  using  only  a  few  layers  of  finite  elements  to  model  the 
1 3 

surrounding  fluid.  Therefore,  for  the  infinite  elements  to  be  competitive  with  finite  elements,  comparable 
results  should  be  obtained  using  infinite  elements  only,  or  much  better  results  should  be  obtained  by  using 
infinite  elements  with  one  layer  of  finite  elements. 

The  computations  for  carrying  out  the  evaluations  of  the  infinite  elements  were  made  using  the  finite 
element  program  NASTRAN.  Separate  programs  generate  the  matrices  for  the  infinite  elements,  and 
NASTRAN  reads  the  matrix  entries  during  execution. 

As  a  preliminary  test  to  validate  the  formulation  of  the  infinite  elements,  the  elements  were  used  to 
compute  the  solution  to  aDinchlet  problem  y2u  —  0  outside  a  boundary  formed  by  a  cylindrical  shell  with 
finite  length  and  closed  ends.  The  Laplace  equation  describes  the  acoustic  field  in  the  added  mass 
approximation.  This  test  evaluates  the  performance  of  the  infinite  elements  for  representing  the  Laplace 
operator  while  avoiding  the  complications  of  the  fluid-structure  interaction  problem  Solutions  obtained  using 
infinite  elements  were  compared  with  solutions  obtained  using  only  finite  elements.  The  standard  for  the 


j\jl R257juuuvift uwiaju'i kVw»uww  vwvu  sr~5r*nmirv.  »v  jr«  jfu  jr^ra^iTJ  jraxu  s  wtwvws  it  ji  w\  u-*t  kjx  xttsa 


comparisons  was  a  solution  obtained  by  the  boundary  element  program  BEASY.  Details  of  the  evaluation  and 
results  of  this  test  are  given  in  Appendix  A.  The  results  of  this  test  show  that,  when  modeling  the  exterior 
region  for  the  three-dimensional  Laplace  equation,  the  use  of  infinite  elements  produces  more  accurate 
solutions  than  the  use  of  finite  elements. 

A  factor  in  determining  the  relative  efficiencies  of  two  methods  is  the  expense  of  setting  up  and  solving 
a  problem  using  each  method.  Since  the  infinite  and  finite  elements  used  in  this  work  each  contain  eight 
nodes,  the  effort  needed  to  produce  a  finite  element  model  and  to  perform  the  computations  using  a  layer  of 
infinite  elements  covering  the  surface  of  the  structure  is  approximately  equal  to  the  effort  needed  for  one  layer 
of  finite  elements.  Thus  the  effort  lequired  in  setting  up  and  solving  a  fluid  model  consisting  of  two  layers  of 
finite  fluid  elements  is  about  the  same  as  the  effort  required  for  a  model  consisting  of  infinite  fluid  elements 
over  a  layer  of  finite  fluid  elements. 

Two  problems  were  used  to  evaluate  the  performance  of  the  infinite  elements  in  applications  to  fluid- 
structure  interaction  problems. 

The  first  pioblem  is  to  compute  some  lower  natural  frequencies  of  a  submerged  flat  plate.  The  plate 
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used  in  this  test  had  the  same  dimensions  and  material  properties  as  one  which  Marcus  used  to  demonstrate 
the  application  of  finite  elements  for  computing  natural  frequencies  of  submerged  plates,  lie  showed  that  only 
a  few  layers  of  finite  fluid  elements  are  sufficient  to  compute  accurately  the  natural  frequencies  of  submerged 
plates,  and  he  obtained  good  accuracy  using  three  layers.  Therefore,  for  this  test,  a  set  of  standard  natural 
frequencies  was  obtained  using  a  conserv ative  model  of  the  surrounding  fluid  consisting  of  five  layers  of  finite 
fluid  elements.  The  accuracy  of  the  natural  frequencies  computed  for  the  test  ;ases  was  determined  by 
comparing  those  frequent  ics  with  the  set  of  standard  natural  frequencies.  Three  models  of  the  surrounding 
fluid  were  used  for  computing  the  natural  frequencies  of  the  plate.  The  first  model,  the  most  desirable  from 
the  standpoint  of  ease  of  modeling  and  economy  in  computation,  represents  the  surrounding  fluid  with  only 
infinite  elements.  The  second  model  uses  infinite  elements  to  lepiescnt  the  infinite  legion  of  fluid  and  a  layer 
of  finite  fluid  elements  undei  the  infinite  elements  to  obtain  better  resolution  of  the  sound  field  in  the  fluid 
adjacent  to  the  structure.  To  get  this  improved  resolution,  one  must  invest  additional  effort  in  preparing  the 
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finite  element  model  and  in  computing  the  solution.  The  third  model  of  the  surrounding  fluid  consists  of  two 
layers  of  finite  fluid  elements  and  no  infinite  elements.  This  model  includes  the  same  number  of  elements  to 
be  constructed  and  (except  for  some  boundary  constraints)  the  same  number  of  degrees  of  freedom  as  the 
second  model.  Thus  the  third  model  requires  nearly  the  same  investment  in  modeling  effort  and 
computational  expense  as  the  second.  Appendix  B  gives  details  of  the  model  used  and  techniques  of  the 
computation. 

Table  1  gives  the  frequencies  computed  in  these  tests. 


TABLE  1.  Natural  frequencies  of  a  submerged  plate. 


Layers 
of  Finite 
Elements 

Infinite 

Elements 

Used 

Frequencies(  normalized) 

i 

2 

Mode 

3 

4 

5 

5 

No 

1.000 

1.000 

1.000 

1.000 

1.000 

0 

Yes 

1.044 

1.024 

1.033 

1.016 

1.020 

1. 

Yes 

1.024 

1.024 

1.021 

1.006 

1.015 

2 

No 

1.052 

1.024 

1.036 

1.027 

1.036 

In  this  table,  the  natural  frequencies  have  been  normalized  by  dividing  the  frequencies  for  each  mode  by  the 
standard  frequency  for  that  mode.  (The  unnormalized  frequencies  are  tabulated  in  Appendix  B)  Table  1 
shows  that  using  only  infinite  elements  to  represent  the  surrounding  fluid  produces  good  results  and  that 
adding  a  layer  of  finite  elements  under  the  infinite  elements  produces  better  results.  However,  it  also  shows 
that  using  infinite  elements  does  not  give  a  significant  advantage  over  using  only  finite  elements,  since  using  a 
layer  of  finite  elements  and  infinite  elements  requires  approximately  the  same  modeling  and  computational 
effort  as  using  two  layers  of  finite  elements. 

The  second  problem  for  evaluating  the  performance  of  the  infinite  elements  is  to  compute  the 
frequency  response  of  a  submerged  cylindrical  shell.  The  cylindrical  shell  had  a  finite  length  and  was 
subjected  to  harmonic  excitation  through  a  range  of  frequencies.  It  was  driven  at  one  point  by  a  sinusoidal 
force  of  amplitude  1.0  and  the  response  was  computed  at  a  separate  transfer  point.  Appendix  B  gives  details 
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of  this  test.  As  in  the  preceding  problem,  the  standard  for  comparison  is  the  solution  obtained  when  the  fluid 
was  modeled  using  five  layers  of  finite  fluid  elements.  Figures  3  to  6  show  comparisons  of  responses  obtained 
using  infinite  fluid  elements  alone  or  over  a  layer  of  finite  fluid  elements  with  the  responses  obtained  using 
only  one  or  two  layers  of  finite  fluid  elements.  The  response  shown  m  these  figures  is  impedance.  The 
impedance  is  given  in  units  of  lb-sec/in  and  equals  the  force  applied  at  the  drive  point  divided  by  the  velocity 
of  the  transfer  point,  the  point  at  which  the  response  is  computed.  Figures  2  and  3  show  that  the  use  of  only 
infinite  fluid  elements  does  not  yield  results  as  accurate  as  does  the  use  of  either  one  or  two  layers  of  finite 
fluid  elements.  Figure  4  shows  that  the  results  obtained  using  infinite  fluid  elements  over  a  layer  of  finite 
fluid  elements  are  accurate  and  are  an  improvement  over  the  results  produced  using  one  layer  of  finite  fluid 
elements.  Figure  5,  however,  shows  that  the  same  configuration  using  infinite  fluid  elements  produces  no 
more  accurate  results  than  does  using  an  equally  costly  model  that  contains  two  layers  of  finite  fluid  elements. 


DISCUSSION  OF  TEST  RESULTS 

The  evaluations  of  the  infinite  elements  show  that  they  produce  accurate  lesults,  especially  if  they  are 
used  over  a  layer  of  finite  fluid  elements.  This  conclusion  is  consistent  with  earlier  work  on  two-dimensional 
infinite  elements  in  which  good  accuracy  was  obtained  when  finite  fluid  elements  were  used  with  infinite 
elements.^  However,  when  solutions  obtained  using  infinite  elements  are  compared  with  solutions  obtained 
using  only  finite  elements,  taking  into  consideration  both  the  effort  to  construct  the  model  and  the  accuracy  of 
the  solutions,  it  is  seen  that  for  added  mass  problems  infinite  elements  do  not  offer  an  advantage  over  finite 
elements.  For  a  more  complete  evaluation  of  the  infinite  elements,  elements  with  higher  order  mapping  and 
shape  functions  could  be  tested.  But  for  a  fair  comparison,  a  higher  order  infinite  element  must  compete  with 
an  equally  higher  order  finite  element  so  that,  as  the  accuracy  of  the  infinite  elements  improves,  so  also  will 
the  accuracy  of  its  competitor.  It  was  judged  that,  for  added  mass  problems,  the  prospects  that  the  accuracy  of 
the  infinite  elements  would  improve  more  than  tnat  of  the  competing  method  was  not  good  enough  to  warrant 


the  investment  of  additional  study,  and  there  are  no  plans  to  continue  the  investigation  of  infinite  elements 


for  added  mass  problems. 


INFINITE  FIND  ONE  LAYER  OF  FINITE  FLUID  ELEMENTS 
FINITE  FLUID  ELEMENTS  ONLY  -  1  LAYER 
FINITE  FLUID  ELEMENTS  ONLY  -  5  LAYERS 


It  is  expected  that  infinite  elements  hold  more  promise  for  applications  to  the  Helmholtz  equation 
problem.  It  is  known  that  a  finite  fluid  element  model  for  the  Helmholtz  equation  case  must  include  six  to 
eight  layers  of  finite  elements  to  follow  the  shape  of  the  wave  and  must  include  damping  boundary  conditions 
on  the  outer  boundary  to  account  for  energy  dissipation  by  outgoing  sound  waves.  Infinite  elements 
developed  for  the  Helmholtz  equation  problem  will  include  wave  characteristics  of  the  sound  field  in  their 
shape  functions  and  are  expected  to  account  for  energy  dissipation.  The  results  obtained  with  the  added  mass 
elements  show  promise  that  the  Helmholtz  equation  elements  will  be  able  to  produce  accurate  solutions  with 
the  use  of  one  layer  of  finite  fluid  elements.  If  this  promise  is  realized,  the  use  of  infinite  fluid  elements  will 
offer  a  considerable  advantage  over  the  use  of  only  finite  fluid  elements.  The  infinite  elements  for  the 
Helmholtz  equation  must  also  compete  with  the  boundary  integral  equation  method  This  method  produces 
accurate  solutions,  but  it  is  also  expensive,  since  it  requires  the  solution  of  systems  of  linear  equations  that  are 
not  banded  and  not  symmetric.  Infinite  fluid  elements  offer  the  possibility  that  for  some  problems,  the 
accuracy  needed  may  be  obtained  at  less  expense. 


APPENDIX  A.  VERIFICATION  PROBLEM 


To  evaluate  the  performance  of  the  infinite  elements  for  representing  the  Laplace  operator,  the 
elements  were  used  to  solve  the  boundary  value  problem  for  the  potential  <j> 

S72<f>  —  0 

in  the  region  cuiside  of  a  bounded  surface  with  boundary  conditions  <j>  =  f  prescribed  on  the  surface.  The 
boundary  of  the  exterior  region  was  a  cylindrical  surface  with  finite  length  and  closed  ends.  The  cylindrical 
boundary  has  length  =  50  in.  and  radius  =  10  in. 

Symmetries  in  the  problem  reduced  the  number  of  elements  needed  to  model  the  cylindrical  boundary 
and  the  exterior  region  A  plane  of  symmetry  midway  between  the  ends  and  perpendicular  to  the  axis  of  the 
cylindrical  boundary  permitted  modeling  only  the  half  of  the  length  of  the  cylindrical  boundary  and  the  half  of 
the  region  on  one  side  of  the  plane.  Two  additional  planes  of  symmetry,  mutually  perpendicular  and  each 
containing  the  axis  of  the  cylindrical  boundary,  permitted  modeling  one-fourth  of  the  remaining  cylindrical 
boundary  and  exterior  region.  Therefore,  only  one-eighth  of  the  cylindrical  boundary  and  exterior  region 
needed  to  be  modeled. 

Figure  7  gi-e-.  did  and  side  views  of  the  cylinder  and  the  surrounding  acoustic  fluid.  Three  varieties  of 
infinite  elements  were  used  to  model  the  acoustic  fluid  in  this  problem.  The  infinite  elements  on  the  sides 
and  on  the  end,  except  the  circle  of  elements  on  'lie  axis  of  the  cylinder,  are  the  eight-node  infinite  elements 
dcscrP  ed  earlier  in  this  report.  In  the  corner  region,  between  the  side  and  end  sections  of  elements,  the 
infinite  elements  have  been  modified  to  extend  to  infinity  in  two  directions,  r  and  z.  The  circle  of  infinite 
elements  on  the  axis  of  the  cylinder  at  the  uid  of  the  cylindiital  shell  are  similar  to  the  eight-node  elements, 
but  their  bases  are  triangles. 

The  infinite  elements  were  evaluated  by  cornpaiing  the  accuracy  of  a  solution  obtained  using  a  model 
that  incorporated  infinite  elements  with  the  acuiiaty  obtained  using  models  that  contained  only  futile 
elements.  The  three  finite  element  models  co  luiined  (wo,  tin ce,  oi  four  layeis  of  finite  elements.  Since  the 
potential  becomes  zero  at  large  radii,  the  bouiulaiy  condition  4>  —  0  was  assigned  at  the  outei  boundary  of  the 
f  .ne  element  model.  For  the  infinite  element  solution,  the  model  of  the  ex'eiio1'  region  consisted  of  one 


Fig.  7.  Finite  and  infinite  element  model  of  the  region  surrounding  the  cylindrical  boundary. 


layer  of  finite  Laplace  elements  over  the  cylindrical  surface  with  infinite  elements  placed  over  this  layer.  This 
evaluation  consists  of  comparisons  of  values  obtained  for  the  exterioi  potential  at  nodes  on  the  suiface 
between  the  finite  and  infinite  elements.  The  standard  for  compaiisons  was  a  solution  obtained  using  the 
boundary  integral  equation  program  BEASY.  The  boundary  conditions  were  <j>(0,z)  —  (l-t/,)cos0  on  the 
cylindrical  surface  and  <f>{0, r)  —  2.6  r  cos 0  on  the  end.  To  be  consistent  with  the  boundary  elements  used 
with  BEASY,  the  boundary  conditions  were  computed  at  0  —  0  ",  15  °,  and  90“  and  linearly  interpolated 
between.  Table  2  gives  potentials  computed  at  several  points,  the  potentials  for  each  point  in  the  table  have 
been  normalized  by  dividing  by  the  standard  potential  computed  at  the  point. 
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TABLE  2.  Potentials  computed  outside  the  cylindrical  shell. 


Point  Coordinates 
r  6  z 

Potentialsf  normalized) 

BIE 

IE 

Methods 

FE4 

FE3 

FE2 

6.67 

22.5 

29.17 

1.000 

0.955 

1.050 

1.044 

1.027 

10.0 

45.0 

29.17 

1.000 

0.914 

0.875 

0.861 

0.820 

13.14 

56.25 

20.83 

1.000 

0.999 

1.002 

0.979 

0.917 

13.14 

56.25 

8.33 

1.000 

0.990 

0.966 

0.927 

0.840 

Note:  In  Table  2,  the  acoustic  fluid  is  modeled  by 

BIE  -  boundary  integral  equation  program  BEASY 
IE  -  infinite  elements  over  one  layer  of  finite  elements 
_ FEi  -  finite  elements  only,  i  layers _ 
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APPENDIX  B.  EVALUATION  PROBLEMS 

Two  test  problems  were  used  to  evaluate  the  use  of  infinite  elements  for  applications  to  fluid-structure 
interaction  problems. 

The  first  test  problem  is  to  compute  the  natural  frequencies  of  a  submerged  cantilevered  plate  This 

13 

plate  has  the  properties  of  the  first  of  two  plates  analyzed  by  Marcus,  and  the  problem  corresponds  to  his 
deeply  submerged  case.  The  length-to-choid  ratio  of  the  plate  equals  2.0  and  the  thickness-to-chord  ratio 
equals  0.0131.  The  plate  has  the  material  properties  of  steel,  with  Young’s  modulus  —  3  0  107  lb/in2, 
Poisson’s  ratio  =  0.3,  and  mass  density  -=  7.324  10' 4  lb-  sec2/in4.  The  fluid  is  assumed  to  have  mass  density 
—  9.34  10  5  lb-  sec2/  in*1  and  sound  speed  —  6000  in/sec.  The  frequencies  computed  here  compare  well  with 
those  computed  by  Marcus.  However,  his  formulation  is  slightly  different  from  the  formulation  used  here, 
since  he  included  the  fluid  matrix  Q  that  derives  from  the  d2p/cft2  term  (see  Eq.  2),  so  his  results  are  not 
compared  with  the  results  from  this  work.  Rather,  comparisons  were  made  using  a  standard  consistent  with 
the  formulation  of  the  added  mass  approximation  used  in  this  work. 

By  using  boundary  conditions  that  represent  the  line  of  symmetry  running  lengthwise  along  the  center 


_  « 


r 


of  the  plate,  o  ’.y  one-half  of  the  plate  in  the  chordwisc  direction  need  be  modeled.  The  finite  clement  model 
of  the  plate  consists  of  five  NASTRAN  QUADl  plate  elements  along  the  length  and  three  elements  across  the 
width.  One  end  of  the  plate  is  constrained  to  represent  \e  clamped  end  condition.  The  fluid  is  modeled  with 
NASTRAN  three-dimensional  IIIEXl  isoparametric  finite  elements  and  with  infinite  elements.  The  fluid 
elements  extend  for  two  elements  outside  of  the  plate  in  the  directions  of  the  length  and  the  width.  using 
half  of  the  mass  and  bending  stiffness  of  the  plate,  only  the  fluid  on  one  side  need  be  explicitly  modeled.  On 
the  outer  boundary  of  the  fluid  finite  element  model  the  boundary  condition  p  =  0  is  applied.  Figure  8  shows 
the  finite  element  model  of  the  plate  and  the  acoustic  fluid. 

The  finite  clement  solution  of  the  submerged  plate  problem  is  implemented  using  the  NASTRAN 
program  in  three  runs.  The  first  NASTRAN  run  computes  a  force  vector  using  a  unit  pressuic  load  on  the 
fluid-structure  interface.  The  resulting  pressure  load  vector  is  processed  to  produce  the  area  matrix  A,  and  in 
the  second  NASTRAN  run  the  matrix  A  is  used  to  foim  the  added  mass  matrix.  The  third  NASTRAN  run 
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Fig.  8.  Finite  clement  model  of  the  submerged  plate  and  surrounding  fluid. 

adds  the  added  mass  matrix  to  the  mass  matrix  of  the  structure  and  computes  the  natural  frequencies  of  the 
fluid-structure  system. 

13 

Several  configurations  of  fluid  elements  were  used  to  compute  the  natural  frequencies  Since  Marcus 
showed  that  for  low  frequencies  the  natural  frequencies  of  a  submerged  plate  can  be  computed  accurately 
using  only  a  few  layers  of  finite  elements  (he  got  good  accuracy  using  three  layers),  a  set  of  natural 
frequencies  was  computed  using  five  layers  of  finite  fluid  elements  and  this  set  was  used  as  the  standard  for 
comparison.  Comparisons  with  frequencies  obtained  using  thiee  and  four  layers  showed  that  the  conservative 
model  with  five  layers  v  as  well  converged.  Table  3  shows  frequencies  obtained  in  the  test. 

'Hie  second  test  problem  is  to  compute  the  forced  frequency  response  of  a  submerged  cylindrical  she'll 
with  closed  ends.  The  amplitude  of  the  response  is  computed  at  a  transfer  point  on  the  cylinder  for  a  range  of 
frequencies.  The  response  is  due  to  a  harmonic  excitation  of  unit  amplitude  at  a. separate  drive  point 
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TABLE  3.  Natural  frequencies  of  a  submerged  plate. 


Frequencies(Hz) 

Layers 
of  Finite 
Elements 

innniUi 

Elements 

Used 

Mode 

1  2  3  4  5 

5 

No 

5.04  28.7  32.9  94.9  95.5 

0 

Yes 

5.26  29.4  34.0  96.4  97.4 

1 

Yes 

5.16  29.4  33.6  95.5  96.9 

2 

No 

5.30  29.4  34.1  97.5  98.9 

The  finite  element,  model  for  the  surrounding  fluid,  in  which  infinite  elements  are  placed  over  one  layer 
of  finite  elements,  had  the  same  configuration  as  that  used  for  the  exterior  region  in  the  Dirichlet  problem  in 
Appendix  A  (see  Fig.  7).  As  in  the  Dirichlet  problem,  symmetries  reduced  the  number  of  elements  needed 
to  model  the  cylinder  and  the  surrounding  fluid,  so  that  only  one-eighth  of  the  cylinder  and  surrounding  fluid 
needed  to  be  modeled.  The  finite  element  analysis  for  this  problem  is  also  implemented  using  the  NASTRAN 
program.  The  matrices  for  the  infinite  elements  and  the  area  matrix  are  generated  by  a  separate  program  and 
entered  during  the  NASTRAN  run  on  DMIG  cards  images.  As  in  the  submerged  plate  problem,  the  structure 
(that  is,  the  cylindrical  shell)  was  modeled  with  NASTRAN  QUADl  plate  elements  and  the  fluid  was  modeled 
with  NASTRAN  three-dimensional  HIEXl  isoparametric  finite  elements.  Figures  3-6  show  frequency 
response  curves  for  the  submerged  cylindrical  shell. 
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